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Abstract 

We study the effect of damping on the generation of baryon asymmetry of the 
Universe in the standard model of the eletroweak theory with simple extensions 
of the Higgs sector. The propagation of quarks of masses up to about 5 GeV are 
considered, taking into account their markedly different dispersion relations due to 
interaction with the hot electroweak plasma. It is argued that the contribution of 
the b quark can be comparable to that of the t quark calculated earlier. 
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1 Introduction 



The discovery of baryon number violation in the standard electroweak theory has led 
to the possibility of constructing a scenario for baryogenesis at the electroweak scale in 
the early universe. This violation, although exceedingly small at the present epoch, can 
be unsuppressed at the then prevailing high temperature The other two Sakharov 
conditions 0] could also be met: C and CP violation originate from the interaction of 
quarks with the Higgs fields. Also, if the electroweak phase transition is of first order, 
the motion of the walls of bubbles of broken phase within the unbroken medium would 
produce the required departure from thermal equlibrium. 

A natural and generic mechanism for electroweak baryogenesis was proposed by Cohen, 
Kaplan and Nelson (CKN) |^. Here a CP-odd charge (like the lepton number) is 
separated by the reflection and transmission of fermions by the bubble wall. It is then 
converted into an asymmetry in baryon number by the baryon number violating processes 
occuring outside the bubble. As a straightforward calculation of the observed baryon 
asymmetry requires too big a CP violation to be available in the minimal standard model 
(MSM), these authors chose to work with simple extensions of the minimal version. The 
required CP violation is achieved by the complex space dependent fermionic mass function 
within the bubble wall, which arise quite generally in such models 0. 

Shaposhnikov [0 studied a similar mechanism of direct separation of baryonic num- 
ber by the bubble wall, taking into account the effect of the hot ambient plasma on the 
quark propagation. As a result of interaction with the quanta of the medium, the quarks 
acquire temperature dependent effective masses and satisfy altered dispersion relations. 
The corresponding fermionic modes or quasiparticles have very different reflection and 
transmission coefficients in different regions of momenta. The details are worked out in 
Farrar and Shaposhnikov (FS) f^. It is found that quark momenta relevant for baryogen- 
esis are much lower than the temperature of phase transition. As a consequence, the CP 
violation of the MSM model suffices to generate the observed baryon asymmetry of the 
Universe. 

Recently Gavela et al and Huet and Sather object to the FS analysis, pointing 
out that they do not include the imaginary part of the effective quark mass giving rise to 
damping of amplitudes. Including the damping they find the reflection coefficients and 
hence the baryon asymmetry to reduce to a negligible value. 



This objection is not agreed upon by FS [O, however. According to them, although 



the leading quark interaction with the medium can be simulated by the real effective mass 
with an alterted dispersion relation, higher order interaction giving rise to the imaginary 
part in the mass and hence damping cannot be treated in this way, as it gives non-unitary 
description. Instead, such dissipative processes must be treated in a many body context. 

Pending a proper discussion of the many body problem, it is useful to study the 
dissipative quantum mechanical problem itself in connection with the baryon asymmetry. 
Here we consider simple extensions of MSM following the work of CKN. However, while 
they considered the t quark propagation, we consider the propagation of lighter quarks, for 
which the finite temperature corrections are important. The reflection and transmission 
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amplitudes are obtained in an iterative series in powers of the mass parameters. The 
results for reflection and transmission currents confirm the loss of unitarity for finite wall 
thickness. However, as the wall thickness tends to zero, unitarity is restored, but the 
fermion-antifermion asymmetry is lost in the model. We reintroduce the asymmetry by 
assuming a very large imaginary part of the mass function within the wall and carry out 
a representative calculation of the baryon asymmetry. 

In sec. 2 we review the propagation properties of quark excitations in the electroweak 
plasma taking damping into account. In sec. 3 we solve the Dirac equation within 
the bubble wall where the mass function is space dependent. It provides the matching 
condition to be used in sec. 4 to find the reflection and transmission amplitudes. In sec. 
5 we calculate the baryon asymmetry in the thin wall approximation but retaining the 
imaginary part of the mass function. Finally in sec. 6 we conclude with a discussion of 
the results obtained. 



2 Quark propagation in hot plasma 

In this section we collect the properties of light quark excitations as they propagate 
through the electroweak plasma at about the phase transition temperature. The most 
important effect of the medium on the quark is obtained by calculating the quark propaga- 
tor at finite temperature 1^. There arises a temperature dependent, chirally invariant 



complex mass with a modified dispersion relation. Neglecting the electroweak contribu- 
tions compared to that due to strong interaction, the leading contribution to the real 
part Eq and imaginary part 7 of the effective mass are the same for both the left(L)- and 
right (i?)-handed quasiparticles, 

Eo = (27ra,/3)^/2T ^ .5T 

and 

7 = .15asT ^ .2T 

with as = .12 at the Z boson mass. For excitations close to Eq, the effective Lagrangian 
incorporating the altered dispersion relation is p|, P, Il0[ 



+ m*R^L, (1) 

3 3 

where we have also included the quark mass acquired through Higgs mechanism. 

Having incorporated the effects of field degrees of freedom in (1), the problem reduces 
to quantum mechanics of left- and right-handed quasiparticles, having the structure of a 
resonance of width 7. We prefer to make the momentum variable complex rather than 
the energy, so that the spatial propagation will be damped. 
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In the following we consider the one dimensional problem where the quasiparticles 
propagate along the z-axis, normal to the bubble wall. Writing 





R = 

the equations of motion derived from the Lagrangian (1) split into two independent sets. 
Defining 

and considering solutions of positive energy E, they are 
where 

QW^sf^-^'t'^ , (3) 

\ -m(z) -(E-E„ + ii)) * ' 

and a similar one for $' with m replaced by its complex conjugate. It suffices for us to 
work with $ only. Note that the current along the z-axis carrried by the components ipi 
and ip^ of $ is 

i, = $^3$, (4) 

o"3 being the third Pauli matrix. 

The planar bubble wall has a finite thickness, extending from z = to z = zq [|T^. It 
separates the broken phase {z > zo) from the unbroken phase {z < 0). The real part of 
the Higgs induced mass m{z) rises from zero in the unbroken phase through the bubble 
wall to the (almost) zero temperature mass mo in the broken phase. The imaginary part 
is non-zero only within the bubble wall. Their actual shapes will be conveniently chosen 
later in sec. 3 below. 

In the unbroken phase (m = 0), the components ipi and ips (also ip2 and '?/'4) decouple. 
Consider (damped) plane waves along z direction, 3 ~ e*^^, K = k + iT^, k > 0. 
They satisfy the dispersion relations jl^ . 



^± = Eo±^,r„ = ±37, (5) 

the (+) and(-) relations holding for ipi and ip^ respectively. In contrast to the situation 
for a free massless particle at zero temperature satisfying E± = ±k, here a part of the (-) 
branch (0 < E_ < Eq) in (5) is also available for quasiparticle propagation with positive 
energy. The (±) branches are called the normal and abnormal ones respectively. Unlike 
the energy E±, the variable k, however, does not represent the true momentum k± of the 
excitation, the latter being given by 

k± = ±^Ei (6) 
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The reversal of sign of k- in the abnormal branch is in conformity with the same for the 
group velocity, 

Thus the above solution for ipi (ips), belonging to the normal (abnormal) branch, advances 
in the positive (negative) 2;-direction. For solutions ipi^s ~ e^*-^^, ipi and ip3 belong to 
opposite branches. Note that the propagation is always damped. 

In the broken phase, the propagation properties are similar except that the Higgs mo 
couples the two components in $ (and in $'). Again consider damped plane wave along 
2;-direction, $ ~ xe*^^, P — p + iTi^p > 0. The spinor x satisfies 

(3{E-Eo + n)-P 3mo ] =o (g) 

[ -3mo -SiE -Eo + i-f)-P)^ ^ ^ 

We get the dispersion relations by setting the determinant of the 2x2 matrix equal to 
zero. Separating the real and imaginary parts, we get 

(E-E„r-(p/zr^^^^, n=?2iE^. (9) 

The normal (+) and abnormal (-) branches arise on taking the square root, 

E = Eo±^g{p), r = ±7^(p). 



where g{p) = -^1 + ^2^p/^^2 ■ Note that the presence of damping (7 7^ 0) removes any gap 
between the normal and the abnormal branches, which exists for 7 = 0, the dispersion 



relation then reducing to E = E()± yp'^/9 + ttIq. The spinor x is obtained by solving (8). 
Normalizing to unit current, we get 

x^CJ. (10) 

where the components c and s stand for 

c = cosh eJ'l', s = sinh eJ'^ (11) 

where 

cosh^ 



E-Eo+p/3 ^4,^ ^ p/3 + ^7 
\ p/3 ' p/3-i7 



The presence of non-zero damping also brings in the phase (j). 

We note here for later use the Lorentz invariant expression for the density of fermionic 
excitations, 

n — {exp Pp • V + 
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where (3 is the inverse temperature of the fluid in the frame where it is at rest, p'^ is the 
energy-momentum 4-vector of the excitation and v^, the 4-velocity of the medium. In 
the wall rest frame, — {E,p),v'^ — "y{l,v) where 7 = l/Vl — and p is the true 
momentum also given by (6). For p along the positive z-direction, we thus have in this 
frame, p-v = E±{l^v/3), up to linear term in v. In the following we require the densities 
of particles with true momenta towards the wall. In the unbroken phase these are given 
by ^ 

for the (±) modes respectively. In the broken phase the corresponding quantities n^. are 
given by the same expressions with the reversal of sign of v. 



3 Solution inside the bubble wall 

We solve equation (2) in a perturbation series in the mass function. Assume 

^z) = e^'^^-^°+''^>^'^{z), 0<z<zo 

^'(2;) will then satisfy 

— = iR{z)^. 
dz 

R{z) has only off-diagonal elements, 

M{z)\ 



R(z)^ 



M{z) 



where M{z) = 3m{z)e^^^^ Eo+n)^^ ^j^d the tilde stands for complex conjugation and 
change in the sign of 7. We now convert it into an integral equation, 

^(z) = $(0) + i r R(z'W(z')dz' 
Jo 

It has an iterative solution, ^(-z) = l^{z)^{0), where 

E(2) = l + i r dz'R{z') - r dz' r dz"R{z')R{z") + ■■■ 
Jo Jo Jo 

We shall actually need the solution for ^(z) at 2; = zq, 

$(2o) = e3^(^-^°+*^)^°'^^E(2o)$(2o) = n{zoMzo). (13) 

Writing 



we get 



rzo rz' 

a = F(l + / dz' / dz"M{z')M{z") + • • •) 



Jo Jo 



(14) 



(5 = iF{ dz'M{z') + ---) 



Jo 



(15) 



with F = e-37^og3i(i7-Eo)^o Note that in the absence of damping (7 = 0), the tilde 
operation reduces to complex conjugation. 

As already mentioned, simple extensions of the Higgs sector of the standard model can 
provide an additional source of large CP violation for baryogenesis In the standard model 
with a single Higgs doublet, the expectation value of the Higgs field is real everywhere 
during the phase transition. But in multi- Higgs models, some of the components acquire 
space dependent values within the bubble wall. This in turn leads to complex space 
dependent mass function for the quarks having Yukawa couplings to those multiplets. 
These are, in principle, calculable from the model considered but, in practice, will depend 
on the many Higgs self-couplings. Here we avoid this problem by assuming a simple but 
anticipated form the mass function. 



within the bubble wall. The parameter 6 relates to the CP violation in the model. 

It suffices to work out the perturbation series to second order to get the leading 
contribution to the asymmetry in the baryonic currents. With the parametrization (16), 
a and /3 in (14-15) can be obtained explicitly to this order as 



m{z) = — z + i—z{zo — z), 



Zo Zq 



(16) 



a = F(l + 9z^C/ + • • •) 
(3 = iF{?>zoV + 



(17) 



U is quadratic in mo and 5 and V is linear. 



U = Ami + iBmoS + CS"^ 
V = arriQ — ibS. 



(18) 



Each of the coefficients A, a etc. depend only on a = 62:0(7 — i{E — Eq)}, 



1 1 \ 1 1 1 



4 Damping in reflection and transmission 



The damping rate has its origin in the scattering of the fermion under consideration by 
other particles in the plasma. Although it does not alter the densities of particles in 
the unbroken or broken phase, a quantum mechanical amplitude gets attenuated as it 
traverses the plasma. In the present context only the attenuation within the bubble wall 
is relevant, where the baryon asymmetry is produced. Thus we prepare our states near 
the wall of the bubble ||15| . 



To study damping effects, consider, for example, the quasiparticle propagation in the 
normal mode. We send a right-handed fermion towards the domain wall from the unbroken 
phase. Noting the reversal of chirality after reflection at the wall, the incident wave (of 
unit current a.t z = 0) and a reflected wave of amplitude r, say, is given by 

$(^) = (I) e'''' + (^) e-*^^ z<0 (20) 



0/ Vr 

On the right (broken phase), we have only the transmitted wave of amplitude t, say. From 
(10) we get 

^z)=t( )e'''^'~'°\z>zo (21) 



—s 

Eq.(13) now serves as the matching condition needed to find the reflection and transmis- 
sion amplitudes 

c \ _ fa P\ fl 
—s) 1/9 a 



t 

giving 



r = -{sa + cl3)/D,t=l/D,D = ca + sl3. (22) 

In the absence of damping (7 = 0), the tilde operation on a and (3 reduces to complex 
conjugation and we immediately obtain from (22), 

|r|2+|t|2 = l (^ = 0), (23) 

expressing the equality of currents in the two phases. However, in the general case (7 7^ 0), 
this does not hold. The reflection and transmission coefficients can be found in the general 
case to the second order in the mass variables using the results for a and j3 obtained in 
(17-19). The expression are rather clumsy and will not be presented here. But for large 
damping, i.e. large values of (672:0), the leading behaviour of these coefficients is simple 
to state, 

badly violating the current conservation relation (23). 

Thus the damping causes the calculational scheme to violate unitarity and results de- 
rived from such a scheme may be questioned. A satisfactory scheme can only be found in a 
many body formulation of the problem which incorporates the true mechanism responsible 
for damping. 
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5 Problem of baryogenesis 



Unitarity can be restored in the present framework only if the thickness of the bubble wall 
tends to zero (thin wall approximation) . In this limit the quasiparticles would not have to 
cross any intervenning medium to go from one phase to the other, which is precisely the 
region where the loss in probability due to the damping matters in the present problem. 
As Zq ^ 0, the functions A,B,C,a and b approach constants and we get 



a = l + 9zl(^ml + ^moS + ^S') + 0{z',) (25) 



P = 3zoi^6+'-mo) + Oiz',) (26) 

Accordingly the matrix Q{zo) —>■ 1 and no asymmetry between fermion and antifermion 
currents would result. 

Nevertheless, we plan to calculate the baryon asymmetry in this model for zero wall 
thickness to examine the effect of damping it. For this purpose we imagine that 5 is very 
large so that 6zq goes to a non-zero constant, say, Aq < 1, even if zq is very small. Then 

a = l + ^ + 0(Ag) (27) 

(3 = ^Ao + 0{Al) (28) 

Note that they are real and independent of 7. 

Since baryon non-conservation through sphaleron processes involves the left-handed 
fermions and antifermions, we are interested in calculating only the left-handed baryonic 
currents in the unbroken phase. 

Consider the propagation of quark excitation in the normal mode. We have already 
calculated for finite wall thickness the reflection and transmission amplitude when right- 
handed fermions are incident on the wall from the unbroken phase. We rewrite them for 
thin wall approximation, 

r = —{sa + cP) / D, ^ ~ D = ca + sP. 

With (25) and (26) the reflection coefficients become, 

where 

2 

= i + hml + ^ + -hpmoAo, h'^ = 4 (7^ + (p/3)' 

The incident flux is the same for particles and antiparticles, viz. , where is given 
by (12). Considering both the particles and antiparticles, the net contribution to the 
reflected left handed baryonic current is 



Hk 1 

--<(i?+-i?+) (29) 
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Here and in the following a bar on a reflection or transmission coeflicient denotes the 
corresponding quantity for the antiparticle. It is obtained by solving the same eqn.(2) 
with m replaced by m*. 

Next consider transmission in the unbroken phase due to incidence on the wall from 
the brooken phase. On the left there is only a transmitted wave of amplitude t', say 

^z) = e-'^', z<0 

On the right we have both the incident wave and the reflected wave of amplitude r', say. 

Again the matching condition (13) gives 

f s + r'c \ _ f a /?\ /O 
\-{c + r's)) ~ \f3 a)[t' 

giving 

r' = - {sa + cP) /D, t' ^ (s^ - c^) /D 

Note that the coefficients fail to satisfy the current conservation across the wall, |r'p + 
|f P 7^ 1. The problem here is that the current in the broken phase as calculated from the 
wave function using (4) is not — 1 + |r'|^ but contains cross terms, 

j = -1 + |r'p + {r'(s*c- sc*) + c.c.} 

= -|1 - r'{s*c - sc*)p + |rf {1 + \s*c - sc*\^} 

= -|c2-sT + k?|c'-5T, 

on using |cp — |sp = 1. Thus with correct normalization, these reflection and transmission 
coefficients coincide, as usual, with those calculated for incidence from the unbroken phase. 
The transmitted left handed baryonic current in the unbroken phase is given by 



dpdE+ b 



27r dp 



n'^T^ - TV) (30) 



where dE^/dp is the group velocity in the broken phase. Adding (29) and (30) we get 
the total baryonic current in the normal mode in the unbroken phase due to reflection 
and transmission as 

7^^^^+"^+^^^+"^+^ ^^^^ 

In a similar way we can work out the baryonic current in the abnormal mode in 
unbroken phase to be 

--{nt-nl){T_-f_) (32) 
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where 

T_ = l/|ac'-/?sf 

and c', s' are given by same expressions as for c and s with changes of sign of {E — Eq) 
and of 0. 

It is now easy to make an order of magnitude estimate of these currents. Formally 
the k- integrals in (31) and (32) extend over to oo and to 3-E'o respectively. But the 
integrands are highly damped at higher values of /c, not just becuse of the presence of 
the density functions; the transmission coefficients for a realistic (i.e. smooth and finite 
width) bubble wall would have fallen exponentially for even lower k values. We set Eq as a 
resonable upper limit for both the integrals. Hopefully the low momentum approximation, 
on which the effective Lagrangian (1) is based, would admit this upper limit. 

As the temperature of the phase transition is about 100 GeV, we may expand the 
density functions in v for small v and approximate the exponential by unity to get 

nl-nlc^ -UvE^ (33) 

D 

Thus the currents (31) and (32) become 

_ Bvm^ fEo dk k p 

Observe the large cancellation in the sum of J+ and J_ giving the total CP-violating 
left-handed baryonic current in the unbroken phase. Evaluating the resulting integral 
approximately we get 

Jcp = J+ + J- 

(35) 

The final step is to obtain the baryonic density in the broken phase from the steady 
state solution to the rate equations in the two phases. CKN|Q find numerical solution 
to the Boltzmann equation. We shall follow FS p|, who solve the diffusion equations for 
small bubble wall velocity to get 

UB = JSpf (36) 

where / is a given function of the diffusion coefficients for quarks and leptons, the wall 
velocity and the sphaleron induced baryon number violation rate. Their estimate for / is 
10~^ < / < 1 in MSM, which should also be valid for its simple extensions. 

Noting the one dimensional entropy density s = 737r/3/3, the baryon to entropy ratio 
is obtained as 

ub/s ~ 1.3 X 10"V^moAo//(2 + A^) (37) 

With mo = 5GeV,P = lO''^GeV'\v = 0.1 and Aq = 1, we get ub/s ~ 3 x 10^7, to 
be compared with the observed value n^/s ~ 5 x 10^^^. 
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6 Conclusion 



We have studied the effect of damping on quasiparticle propagation in the electroweak 
plasma in connection with the problem of baryogenesis in non- minimal versions of the 
standard model, where the Higgs sector is extended to include more than one multiplet. 
In general, such an extension gives rise to a complex mass function for a quark within 
the bubble wall formed during the phase transition. This constitutes the CP violation 
needed for baryogenesis. Its real and imaginary parts are parametrized in a simple way 
making the integrals simple to evaluate. We follow closely the technique of CKN0, but 
consider a direct separation of baryon number by the bubble wall rather than of some 
other CP-odd charge 0. 

The inclusion of the temperature dependent effective mass gives rise to two modes of 
quasiparticle propagation in the plasma. Our calculation shows that both modes must 
be taken into account. In fact, the net baryon asymmetry current results after large 
cancellation between the baryonic currents carried separately by the two modes. 

We also include the recently discussed damping suffered by the quasiparticles while 
propagating through the plasma |10|. Like the effective mass, it also arises from in- 
teraction of a fermion with other particles in the plasma. This damping appears to ruin 
an otherwise successful calculation of baryon asymmetry in the MSM. Besides spatial 
attenuation, it also affects the spinor components of $ and the dispersion relations are 
no longer separated by a mass gap. Also the two spinor components acquire equal and 
opposite phase. 

The quasiparticle propagation is damped everywhere within the medium, in the broken 
and unbroken phases, as well as within the bubble wall. However, it is the damping 
associated with the propagation within the wall, which is of concern in the problem of 
baryogenesis. It causes the magnitudes of the reflection coefficients to be reduced, more so 
the larger the value of jzq. The problem one faces here is that the calculational scheme is 
non-unitary and the question of reliability of a calculation within such a scheme remains 
open. 

Unitarity is restored in this formalism only in the thin wall approximation, when a 
calculation of baryon asymmetry would be free from the above objection. But ^1{zq) — * 1 
as ~^ and there is no asymmetry between the fermion and the antifermion as they 
interact with the wall. We are then led to consider a rather unrealistic situation where 
the coefficient S in the imaginary part of the mass function is so large that 6zo = Aq, a 
non- zero constant, even in the thin wall approximation. It is now simple to calculate the 
baryon to entropy ratio of the Universe. As expected, it does not suffer any reduction 
due to damping rate when compared with a similar expression calculated without taking 
it into account 



We now come back to comment on the real problem with finite wall thickness. Though 
a satisfactory framework for calculation should avoid any non-unitarity due to damping 
by including the many body processes responsible for it, it, nevertheless, appears that 
the damping effect in baryon asymmetry will persist even in such a framework. Thus 
the result calculated without damping [O is likely to be reduced by suppression factor 
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appearing in the reflection and transmission coefficients (24). Note that this suppression 
is much less than that encountered in MSM where one has to go to higher order in 
perturbation expansion in the mass matrix. Furthermore, the magnitude of the wall 
thickness in extended versions of the standard model is not known. In view of the many 
parameters in the potential of such a model, the indication is that it could be small in 
several regions of the parameters. For zq < j^GeV~^, the damping factor is of order 10~^ 
at most. Within the uncertainties of the transport properties of the electroweak plasma, 
the baryon to entropy ratio would still be around the observed magnitude. 

CKN H find that the t quark propagation through the bubble wall can produce the 
observed baryon to entropy ratio. We here argue that the contribution of the b quark, 
for which finite temperature correction in mass, in both the real and imaginary parts, is 
important, is also likely to be of similar magnitude. 
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